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ABSTRACT 

This  thesis  models  induced  magnetic  fields  from  the  motion  of  seawater  in  the  earth’s 
magnetic  field  analytically  and  compares  the  results  to  arctic  on-the-ice  magnetic  fluctuation 
measurements.  The  oceans  have  various  types  of  internal  motions,  such  as  internal  waves  and 
turbulence.  This  motion  of  seawater,  which  is  a  conductor,  in  the  earth’s  magnetic  field  induces 
a  current  density.  This  current  density,  in  turn,  induces  its  own  magnetic  field.  This  thesis 
models  internal  waves  and  upper  layer  ocean  turbulence  analytically.  The  corresponding 
induced  magnetic  fields  are  calculated  using  a  static  form  of  Maxwell’s  equations  and  parameters 
for  the  Arctic  are  inserted.  Comparisons  are  made  with  measurements  from  the  A:ctic  Internal 
Wave  Experiment  (1985).  The  predicted  fields  from  internal  waves  have  magnitudes  that  are 
measurable  and  of  the  same  order  of  magnitude  as  ionosphericaliy  generated  fields.  The 
predicted  fields  from  turbulence  are  several  orders  of  magnitude  smaller  than  ionosphericaliy 
generated  fields.  Besides  giving  information  about  internal  motions  in  the  ocean,  the  seawater 
induced  fields  are  a  noise  source  in  magnetic  anomaly  detection. 
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I .  INTRODUCTION 


Seawater  is  a  conductor  situated  in  the  earth's  magnetic 
field.  Motion  of  this  conductor  from  surface  waves,  internal 
waves,  turbulence,  etc.  produces  a  motional  electromotive 
force.  The  resulting  current  density  induces  its  own  magnetic 
field.  Understanding  the  production  and  propagation  of  such 
magnetic  signatures  provides  information  about  the  underlying 
seawater  motion.  Also,  these  induced  magnetic  fields  are  a 
noise  source  for  other  magnetic  measurements,  such  as  magnetic 
anomaly  detection. 

This  thesis  constructed  analytical  models  of  these  induced 
fields  and  compared  the  predicted  values  with  experimental 
measurements.  One  model  considered  internal  wave  sources  and 
another  ocean  upper  boundary  layer  turbulence.  The  internal 
wave  model  had  versions  for  the  Arctic  and  lower  latitude 
oceans. 

The  internal  wave  model  started  with  the  mathematical 
derivation  of  a  governing  equation  for  seawater  velocity  from 
fluid  mechanics.  The  model  assumed  no  rotation,  zero 
viscosity,  incompressibility,  and  slow  variation  of  mean 
density  with  depth.  Separation  of  variables  provided  a  linear 
differential  equation  for  the  amplitude  of  the  vertical 
component  of  velocity.  The  corresponding  solution  and  the 
separation  of  variables  equations  together  provided  an 
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analytical  model  of  seawater  velocity,  v.  The  velocity 
frequency  dependence  was  then  tailored  to  fit  observed 
velocity  spectra  for  either  the  Arctic  Ocean  or  lower  latitude 
oceans . 

Next,  the  induced  magnetic  field  was  calculated.  The 
motion  of  the  conducting  seawater  with  velocity  v  in  the 
earth's  magnetic  field  B  (assumed  constant)  induces  an 
electromotive  force  and  causes  a  current  density  J.  This 
current  density  J  induces  its  own  magnetic  field  B' .  Because 
the  frequencies  involved  are  on  the  order  of  10'3  H2  or  less, 
displacement  currents  were  neglected  and  the  inductive  field 
B'  calculated  with  the  Biot-Savart  Law.  The  total  induced 
field  B*  at  a  field  point  was  calculated  assuming  an  internal 
wave  of  infinite  y-extent  with  wavenumber  in  the  x-z  plane. 
Integration  over  all  x  and  y,  and  over  z  within  a  depth  range 
from  D0  to  D  yielded  B' .  The  wavenumber  dependence  was 
removed  by  two  different  methods.  In  one,  B'  was  converted  to 
a  power  spectrum,  transformed  from  a  one-dimensional  to  a  two- 
dimensional  spectrum,  and  integrated  over  wavenumber  from  a 
minimum,  kmjn,  to  infinity  to  give  B'2(<.>).  In  the  other, 
application  of  a  boundary  condition  to  the  vertical  component 
of  velocity  produced  a  dispersion  relation  between  o  and  kx. 
Choosing  the  lowest  mode  of  oscillation,  using  the  dispersion 
relation,  and  converting  to  a  power  spectrum  yielded  B,2(o). 

Data  from  the  Arctic  Internal  Wave  Experiment  (AIWEX)  in 
1985  provided  a  means  to  compare  the  model  B'2  power  spectral 
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density  as  a  function  of  frequency  with  experimental  results. 
The  predicted  magnitudes  for  the  induced  fields  B'2  were  in 
the  same  range  as  extrapolated  magnetometer  measurements.  The 
large  spatial  extent  of  internal  waves  produced  signals  of  the 
same  order  of  magnitude  as  ionospherically  generated  signals. 
The  first  k-dependence  method  yielded  frequency  responses  that 
do  not  follow  the  1/f2  dependence  of  the  data.  The  second 
method  does  yielded  the  1/f  or  1/f2  frequency  dependence 
expected.  The  success  of  the  second  method  implied  that  the 
modal  structure  of  the  internal  wave  field  must  always  be 
utilized  in  calculating  the  induced  magnetic  fields. 

The  turbulence  induced  magnetic  field  model  had  a  similar 
formulation.  It  assumed  the  region  of  turbulent  seawater 
resides  in  the  upper  boundary  layer  of  the  ocean  between 
depths  D0  and  D,  and  to  be  three  dimensionally  isotropic 
within  this  region.  The  seawater  velocity  v  was  analytically 
modeled  as  a  linear  superposition  of  plane  waves  propagating 
in  three  dimensions  with  various  wavenumbers  and  frequencies. 
This  motion  in  the  earth's  magnetic  field  B  resulted  in  the 
current  density  J,  which  was  utilized  in  the  Biot-Savart  Law 
to  calculate  the  induced  magnetic  field  B'  at  a  field  point, 
again  neglecting  displacement  current.  The  total  field  was 
found  by  integrating  over  all  x  and  y,  and  over  z  between  D0 
and  D.  Keeping  the  principal  term  of  the  B*  expression,  the 
induced  field  components  were  integrated  over  wavenumber  from 
a  minimum,  k  •  ,  to  infinity.  Then  the  horizontal  field 
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magnitude,  B'H,  was  converted  to  a  power  spectral  density  as 
a  function  of  frequency. 

This  power  spectrum  was  compared  to  the  AIWEX  magnetometer 
data.  The  frequency  dependence  of  the  model  agreed  with  the 
observed  1/f2  dependence  of  the  data.  However,  the  predicted 
magnitudes  were  several  orders  of  magnitude  below  the  data. 
The  AIWEX  signal  in  the  lO^-lO"2  Hz  range  appeared  to  arise 
from  fields  generated  in  the  ionosphere. 


4 


IX.  BACKGROUND 

A.  OCEAN  ENVIRONMENT 

The  conductivity  of  seawater  principally  derives  from 
dissolved  salts.  In  general,  conductivity  decreases  with 
depth.  An  average  value  for  ocean  conductivity  is  4  (S/m) . 

Seawater  density  varies  in  the  ocean.  Temperature  and 
salinity  are  the  main  governing  factors.  Density  as  a 
function  of  depth  generally  displays  three  zones  (Gross, 
1971) .  In  the  surface  zone,  from  the  surface  to  50-100  m,  the 
density  remains  approximately  constant  with  depth  due  to 
mixing  of  the  upper  layer  of  ocean  by  wave  action.  The  next 
zone,  from  50-100  m  to  around  500-1000  m,  forms  the  pycnocline 
zone.  Density  increases  rapidly  with  depth.  Generally,  the 
density  follows  a  monotonic  increase  because  the  temperature 
declines  with  depth.  The  bottom  zone,  below  the  pycnocline, 
is  called  the  deep  zone.  Here,  density  slowly  increases  with 
depth  due  to  decreasing  temperature.  Fig.  l  (Gross,  1971) 
depicts  these  zones  and  their  general  temperature,  salinity, 
and  density  variations. 
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Figure  1:  Typical  Ocean  Thermocline,  Halocline,  and 
Pycnocline  (Gross,  1972) 

The  earth's  permanent  magnetic  field  permeates  the  oceans. 
Since  the  field  is  approximately  that  of  a  magnetic  dipole, 
the  direction  and  magnitude  of  the  field  vary  across  the 
oceans.  It  goes  from  basically  horizontal  with  respect  to  the 
surface  of  the  earth  and  magnitude  of  about  30,000  gamma 
(nanotesla)  near  the  equator,  to  basically  vertical  with 
magnitude  of  about  60,000  gamma  at  the  north  and  south 
magnetic  poles.  Seawater  motions  of  the  ocean  within  this 
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field  induce  the  electromotive  force  and  the  resulting  current 
density  J. 

B.  INTERNAL  WAVES 

A  change  in  density  with  depth  (  dp/dz  <  0  )  allows 
buoyancy  oscillations  within  the  ocean  known  as  internal 
waves.  The  momentum  and  continuity  equations  for  fluids  are 
used  to  derive  an  equation  for  the  vertical  velocity  of  these 
internal  waves. 

Following  Gill  (1982),  the  fluid  momentum  equation  is 

-^+2flxu=-p-lVp-^  +  v^u,  (1) 

where  u  =  fluid  velocity  =  ui  +  vj  +  wk,  fl  «*  system  rotation 
angular  velocity,  and  v  =  kinematic  viscosity.  Assume 
isentropic  motion,  i.e.  no  viscous  effects  (v=0)  and  no 
rotation  (fl=0) .  The  mass  equation  of  continuity  is 

+  V'U  =  0 . 

r  Dt 

Assume  density  is  a  function  of  potential  temperature,  0,  and 
salinity,  s,  independent,  of  pressure,  p.  Using  the  chain  rule 
and  taking  e  and  s  constant  with  respect  to  time  over  a 
differential  element 

Dp  _  +  _d£  _Ds  _  0  (3) 

Dt  ~  M  Dt  ds  Dt 
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The  continuity  equation  then  becomes  the  incompressible  fluid 
condition 


Vu  =  0. 


(4) 


Assume  small  perturbations  in  pressure,  p',  and  density,  p', 

(5) 


P  =  Po  +  p'>  P  =  Po  +  p'f 


dp, 

dz 


£  =  ~gp 


o‘ 


The  fluid  momentum  equation  becomes 


p-ff  -  -vp'  - 


<«) 


Applying  the  convective  derivative  with  average  velocity  v, 
the  time  derivative  of  density  becomes 


M  .  |E  *  <vV)v-  4|^ 
DC  dt  dt 


+  w- 


dp  o 
dz 


=  0 


(7) 


The  fluid  velocity  components  u  and  v  can  be  eliminated  by 
first  taking  the  partial  derivative  with  respect  to  time  of 
the  continuity  equation 


(fu  +  d2  v  +  d2  w  _  q 
dtdx  dtdy  dtdz 


(8) 


and  then  using  Equation  6  in  component  form  to  replace  time 
derivatives  of  u,v  with  spatial  derivatives  of  p*.  The 
result  is 
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(9) 


&W  _  +  &2l_ 

Podzdt  -  dx3  dy2 
Combining  Equation  [7]  and  the  time  derivative  of  the  z- 
coroponent  of  Equation  [6]  yields 


at2 


*  N2w  = 


_i__av 

p0*azat ' 


(10) 


where  N  is  the  Brunt -Vaisala  frequency,  given  by 


N2 


g  dPo 

Po  dz 


(11) 


Taking  the  horizontal  Laplacian  of  Equation  [10] 


dc2  dx 2  dy2  dx2  dy2  Pa  dx2  dy2 


d2p\!2) 

dzdt  ' 


Substituting  Equation  [9]  in  the  right  hand  side  and 
rearranging  yields  an  equation  for  w,  the  z-component  of 
velocity. 


JL{JL. 

at2  dx2 


d2  Id,  a 

-  + - *-  (p  n-T-  )  J  V  + 

dy 2  P  0  dz  dz 


N2(-iL  +  -*L 

dx2  dy 2 


)  w  =  0. (13) 


z, 


Assume  that  the  vertical  velocity,  w,  varies  with  depth, 
much  more  rapidly  than  the  density, p0.  Then 


-4-(  Po4-)w 

p  .  dz  °  dz 


dzw 
dz2  ' 


(14) 
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which  is  known  as  the  Boussinesq  approximation. 
Equation  [13]  then  simplifies  to 


dt2  dx2  dy 2  dz 2  dx2  dy2 


(15) 


This  equation  describes  the  behavior  of  the  vertical  velocity 
component  of  an  internal  wave  disturbance. 

Separation  of  variables  can  simplify  Equation  [15] 

w(x,  y,  z,  t)  =  X(x)  Y(y)  Z(z)  Tit).  d«) 


The  X,  Y,  and  T  portions  have  the  form 

+  kX  =  0,  (17) 

dx2 

with  the  following  solutions 

X  =  X0eik**-,  Y  *  Y0eik'y;  r=roeiwt.  (1*> 

The  z  part  becomes 

&A  +  k2t(l-  Nizy  )  2=0.  (19) 

dz2  w2 

This  equation  possesses  the  Sturm-Liouville  form.  Thus,  its 
solutions  are  orthogonal  eigenfunctions  Z(z)  with  real 
eigenvalues.  The  vertical  velocity  solutions  then  take  the 
form 
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w(x,y,z,  t)  =  Z{z)  ellk**  *  kyV  ~  wC) .  (2°) 

Since  they  are  orthogonal,  any  internal  disturbance  can  be 
expanded  as  a  sum  of  these  solutions. 

For  Z(z) ,  try  solutions  of  the  form 

Z(z)  =  i t(z)  e 1  k *  *.  <21> 


Substituting  this  into  Equation  [19]  gives  for  G( z) 


az 


dz: 


y\)  v  »  o. 


(22) 


The  term  with  N(z)J  makes  this,  in  general,  a  nonlinear 
differential  equation. 

To  simplify  this  equation,  approximate  p(i.)  as  a  piecewise 
continuous  function  over  adjacent  regions  such  that  dp/d z  is 
a  constant  in  each  region.  Then  N(z)  will  be  a  constant  in 
each  region.  Equation  [22]  can  be  solved  for  each  region,  and 
the  solutions  joined  with  boundary  conditions  specifying 
continuity  of  velocity.  Alternately,  approximate  p(z)  by  a 
superposition  of  linear  slope  segments  each  with  dp/dz  equal 
to  a  constant.  Then  N(z)  is  a  constant  for  each  segment  and 
Equation  [22]  is  again  readily  solved.  These  approximations 
are  seldom  applicable  to  the  real  ocean  (Gill,  1982)  for 
accurate  velocity  profiles.  However,  the  resulting  solutions 
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do  have  the  general  shape  and  the  oscillating  character  of  the 
profiles  resulting  from  actual  p(z)  distributions. 

With  N ( z)  equal  to  a  constant,  Equation  [22]  becomes  a 
second  order,  constant  coefficient,  linear  differential 
equation 

+  ptf  =  0,  (23) 

dz2  dz 


where  a  -  2  ikz.  p  ■  [  ( —  )  ( k2x  +  k})  -  Jcf]  .  (24) 


The  solutions  of  the  characteristic  equation  are 


k  ■  .  -i*. .  i 


(  *1.1 Sl ;  ikl  *  ki)  (25) 

U2  y 


Define  6  as 


N 


<  .y  -.-a2  >  ♦  k}) 

CO2 


(26) 


In  a  stratified  ocean,  the  maximum  oscillation  frequency 
equals  the  buoyancy  frequency  N.  Thus,  o  £  N  ,  and  6  is  real. 
For  e  <  N,  the  general  solution  of  equation  [23]  is 

£-(z)  =  e"iJf**[Ci  e<ia*  ♦  C2  e-i4r)  .  (27) 
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Initial  and  boundary  conditions  determine  the  constants  Cl  and 
C2. 

The  vertical  velocity  component  is  then 

#90) 

w(x  ,y,z,t)  =  tf(2)  ei(JCjr*  *  *yy  *  *»*  *  we)  =  v(Z) 

The  incompressible  continuity  condition,  Equation  [6], 
determines  the  horizontal  components,  u  and  v, 

u(xiy,2,c)  =  [i#  -  kt  v)  (  --*■*-  - )  ei(k' T  '  ut) ,  (29) 

**  X 

vU.y.z,  t)  =  [i-§£  -  *3  (  -  -  )  eiik ' r  ‘  uc) .  (30) 

52  +  iCy 

These  analytical  expressions  for  the  seawater  velocity  now 
facilitate  calculating  the  induced  magnetic  fields. 

C.  ARCTIC  UNDER  ICE  SEAWATER  TURBULENCE 

In  the  arctic,  turbulent  seawater  flow  exists  from  below 
the  ice  (effectively  zero  depth)  down  to  some  depth,  D,  at 
which  the  motion  transitions  to  internal  waves.  40  m  is  a 
typical  value  for  D  (private  communication  with  T.  Stanton) . 

Utilize  superposition  of  three-dimensional  plane  waves  of 
varying  wavenumber  and  angular  frequency,  o,  to  model  this 
turbulent  region 
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v-  lvx(kx,o)ik'*£  +  vy(ky.o)elk’yJ  +  vM(k„  u)  eik>'Jc  )  (31) 


Assume  three-dimensional  isotropy  of  turbulence  to  simplify 
the  analysis 

vx(kx,(i>)  -  vy(ky,u)  »  vt(kt,<j))  ■  v(k,«)  .  (32) 


Empirically  model  the  velocity  amplitude  function  v(k,o)  as 


v(k,  o) 


yjkl  *  k2  y/(i>20  +  <i>2 


(33) 


The  power  spectrum  for  this  function  is  then 


v2(ic,oi)  = 


(kl  +  k2)  (u20  +  o2) 


(34) 


The  constant  factors  k0  and  w0  in  the  denominator  keep  the 
power  finite  as  k  and  o  approach  zero.  For  large  k,  this 
exhibits  the  same  k'2  dependence  as  for  internal  waves 
(Garret,  Munk,  1972) .  For  large  u,  this  exhibits  the  o'2 
dependence  characteristic  of  under  ice  internal  wave  fields 
(ONR  Report,  1991) .  The  velocity  expression  is 


v(k,  o),  r,  t)  - 


y/ko  ♦  k2J(D20  *  (j)2 


,Hk  ■  r  -  uc) 


(35) 
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III.  MAGNETIC  SIGNATURES 


A.  PROPAGATION  OF  INDUCED  ELECTROMAGNETIC  FIELDS 

This  model  neglected  attenuation  due  to  the  propagation  of 
the  electromagnetic  waves  through  seawater,  a  conductor, 
because  of  the  low  frequencies  involved  and  the  corresponding 
long  wavelengths.  The  electromagnetic  waves  are  characterized 
by  their  wavelength,  X.  For  a  "good  conductor",  like 
seawater,  and  a  "good  insulator",  like  air,  the  wavelengths 
are 


X 


conductor 


2k 


X 


insulator 


C 

1' 


(36) 


where 

p  =471 X10-7— ,  a  =  cond-  4  —  ,  c=3.0xlD8— .  (37) 

°  A2  in  S 

Table  I  lists  wavelength  values  for  typical  frequencies  of 
interest  in  internal  ocean  seawater  motions. 
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Table  I:  Wavelength  Values 


Fraq  (Ha) 

^•aawatar  * 

*alr  <®>  * 

10° 

1570 

3x10® 

10'1 

5030 

3xl09 

10‘2 

15,700 

3xl010 

10*3 

50,300 

3X1011 

10*4 

157,000 

3Xl012 

The  induced  magnetic  fields  will  "propagate"  with 


characteristic  wavelength  *se#M4tef  up  through  the  ocean  and  then 
with  the  much,  much  longer  wavelength  X8tr  above  the  surface. 

Characteristic  ocean  buoyancy  frequencies  are  around  10'3  - 
10'4  Hz.  Since  internal  wave  frequencies  are  always  less  than 
the  buoyancy  frequency,  and  the  depth  of  the  Arctic  Ocean  is 
around  3000  m  (Dietrich,  et  al.,  1980),  the  corresponding 
wavelengths  of  the  induced  electromagnetic  radiation  are  much 
larger  than  the  depth.  Within  the  ocean  and  in  the  air  above 
the  ocean  surface,  conduction  current,  J,  and  not  displacement 
current,  edE/dt,  creates  the  induced  magnetic  field.  The 
field  point  is  in  the  "near  field"  and  propagation  effects  are 
negligible.  Relevant  turbulence  frequencies  range  around  10’1- 
10'3  Hz.  We  confined  our  interest  to  turbulence  in  the  upper 
boundary  layer  of  the  Arctic  Ocean.  Since  a  nominal  depth  for 
this  layer  is  40  m,  the  electromagnetic  radiation  wavelengths 
involved  are  again  much  larger  than  the  depth  and  fields  above 
the  surface  are  "near  field"  with  negligible  propagation  loss. 
Thus,  we  neglected  attenuation  due  to  propagation  in  the 
seawater  conducting  medium.  For  these  same  reasons,  the 
sea/ice/air  interfaces  do  not  attenuate  the  induced  fields. 


16 


The  earth's  magnetic  field  inside  the  ocean  will  possess 
the  same  value  as  in  the  air  above  the  ocean.  This  follows 
from  the  assumption  that  the  earth's  field  is  a  static  (DC) 
field,  and  that  the  magnetic  susceptibilities  of  the  air  and 
ocean  are  the  same  (i.e.  both  are  non-magnetic) .  The 
equivalence  derives  from  the  resulting  boundary  conditions 
between  the  two  media 

^sea,  tangential  ~  ^air,  tangential  •  ^aea, normal  ~  ^air, normal  ’ 

B.  INTERNAL  WAVE  SIGNATURE 

1.  Geometry 

The  method  of  Crews  and  Futterman  (1962)  allowed  us  to 
calculate  magnetic  signatures.  The  current  density  followed 
from  the  seawater  velocity  in  the  earth's  magnetic  field. 
This  current  density  was  then  utilized  in  the  Biot-Savart 
method  of  calculating  magnetic  induction  at  a  field  point. 

Fig.  2  shows  the  geometry  of  the  earth's  magnetic  field 
within  the  ocean  (northern  hemisphere)  .  To  simplify  the 
calculations,  choose  the  wave  vector  in  the  x-z  plane.  Then 
the  y-components  of  the  wave  vector  and  the  velocity  are  zero. 
Assume  the  wave  extends  infinitely  in  the  ±y  directions.  The 
seawater  in  the  volume  element  dxdydz  at  the  point  (x,y,z) 
with  velocity  v(x,y,z,t)  forms  a  source  element  for  an  induced 
magnetic  field  at  the  field  point  located  at  (0,y,-h).  Fig. 
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3  shows  this  geometry,  "z"  represents  the  depth  below  the  sea 
surface.  "h"  represents  the  height  of  the  field  point  above 
the  sea  surface. 


♦ 

z 


X 


Figure  2:  Geometry  of  the  Earth's  Magnetic  Field 
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Figure  3:  Geometry  for  Calculation  of  Internal  Wave  Induced 
Magnetic  Field 
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2.  Induced  Magnetic  Field 

Seawater  behaves  as  a  relatively  good  conductor,  with 
conductivity  of  approximately  4  S/m.  The  motion  of  this 
conductor  in  the  earth's  magnetic  field  creates  an 
electromotive  force.  This  in  turn  leads  to  a  current  density 
given  by 

o  (  v^B)  ■  o  (  vtBx  -  v^g)  f.  (3F) 


Since  velocity  is  confined  to  the  x-z  plane,  the  current 
density  J  is  in  the  y-direction. 

This  current  density  creates  a  secondary  induced  magnetic 
field.  The  Biot-Savart  law  gives  the  differential  induced 
magnetic  field  due  to  a  current  element  di 

dEi  -  ,  di*JdA.  (40) 

2  n  r2 


The  prime  indicates  the  secondary  induced  field.  Integration 
of  Equation  [40]  gives  the  components  of  the  total  induced 
field.  The  assumption  of  infinite  extent  of  the  wave  in  the 
±y  directions  (i.e.  symmetry  about  y»0)  results  in  the  y- 
component  summing  to  zero.  The  x  and  z  components  are 

sW/dBWX-n-F(ft*z)-  <41> 


(42) 
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Substituting  in  the  expressions  for  di,  v,  B,  and  R  and 
integrating  [see  Appendix  A]  gives 


ij— jj-  )  co»yco»*  ♦ 


k,*l& 


)  ) 


0 


<43) 


i  tt,agc>. 


V.-tat 


[  ( 


^t77T  *  -j^7^r•,co•Yeo,♦  *  *;,in*( 


(U) 

*~ir  ,,e-* 


Here,  6  reduces  to 


W2  - 


(45) 


These  equations  are  fairly  complicated.  However,  they 
simplify  somewhat  when  evaluated  specifically  for  conditions 
in  the  Arctic 


B  8  fl  Jc;  0  *  90°;  cos0  =  0;  sin0  «  1;  Da  »  0. 


(46) 


The  induced  magnetic  field  components  become 


Em?  I  (_1  - is*!—  1  .(47) 

w*  ^78  *,(y'-2wa) 


Simplifying, 


~  )-p-  (  *•«  -  ••a'VJc„cos6£>  -  fisinftD)  )  .  (  49  ) 
Jcj*6a  ** 


Inserting  the  expression  for  6,  these  become 
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Jc,  N*  N  *»*  >|  «*  )|  w* 


(50) 


Note  that  the  x  and  z  components  are  exactly  the  same  except 
for  the  phase  factor  of  i 


B't*  i  B'A 


with  i 


(51) 


Thus,  the  x  and  z  components  are  90*  out  of  phase,  and  the 
induced  magnetic  field  vector  rotates  in  the  x-z  plane.  It 
has  magnitude 

l*'l  =  \b'x\  =  \b‘z\  .  (52) 

Important  characteristics  of  B*  are  (1)  exponential  decrease 
in  magnitude  with  height,  h,  above  the  internal  wave  field, 
and  (2)  increase  in  magnitude  with  the  vertical  height  D  of 
the  internal  wave  field. 

These  expressions  contain  angular  frequency  u,  wavenumber 
kx,  and  time  t.  Comparison  with  experimental  data  often 
requires  the  total  field  power  spectrum  as  a  function  of 
frequency.  Thus,  all  of  the  wavenumber  and  time  dependencies 
must  be  removed. 

The  frequency  and  wavenumber  dependencies  of  the  velocity 
were  empirically  modeled  based  on  actual  measurements  of 
velocity  power  spectra.  The  amplitude  of  the  internal  wave 
velocity,  C,  which  so  far  has  only  been  a  multiplicative 
component  in  the  equations,  became  a  function  of  frequency  and 
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wavenumber.  This,  in  turn,  added  additional  frequency  and 
wavenumber  dependence  to  the  induced  magnetic  field  B'. 

The  Arctic  Internal  Wave  Experiment  (AIWEX) ,  conducted  in 
1985  in  the  Canada  Basin  of  the  Arctic  Ocean  (depth  >  3500  m) , 
measured  horizontal  velocity  with  an  electromagnetic  current 
meter  at  a  depth  of  100  m.  Fig.  4  shows  one  such  spectrum. 
Sampling  rate  limited  the  highest  frequencies  measured  to 
around  6  cycles  per  hour  (1.7xl0*3  Hz),  the  local  buoyancy 
frequency.  The  spectra  indicate  a  1/f  frequency  dependence 
over  the  range  10'1  -  10  cycles  per  hour  (2.8xl0's  -  2.8xl0'3 
Hz)  (Levine,  et  al.,1987).  This  frequency  dependence  was 
empirically  modeled  by 


C2  C2 

|vj2  =  --  ---- .  =  - -  ,  ("A"  for  AIWEX)  .  (53) 

Jo2,  +  u2  2  njfl  +  f2 


For  o  >>  <i>0,  it  behaves  as  l/u.  For  u  <  <*>0,  its  value  stays 
approximately  constant  and  limits  the  total  energy  in  the 
spectrum  as  o  approaches  zero. 

The  Coordinated  Eastern  Arctic  Experiment  (CEAREX)  in  1989 
measured  velocity  power  spectra  at  depths  of  50,  100,150,200, 
and  250  m  below  the  surface  of  the  ice.  The  50  m  measurements 
were  made  with  an  acoustic  doppler  profiler  and  the  other 
measurements  were  made  with  electromagnetic  current  meters. 
Appendix  C  contains  reproductions  of  the  resulting  velocity 
power  spectra.  All  plots  showed  an  approximately  1/f2 
frequency  dependence  in  the  range  0.1-10  cycles  per  hour 
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Figure  4:  AIWEX  Horizontal  Velocity  Spectra  (solid  lines); 
Garrett-Munk  Model  Spectrum  (dashed  lines)  for  comparison. 
(Levine,  et  al. ,  1987) 
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(2.78X10'5  -  2.78X10*3  Hz).  Often  this  dependence  flattened 
out  for  frequencies  below  about  4xl0"2  cycles  per  hour 
(l.llxlO*6  Hz).  (Czipott  and  Podney,  1989). 


This  frequency  dependence  was  empirically  modeled  by 


ve\2  « 


cl 


cl 


o  <,  +  o* 


4n2  (fj  ^  fJ) 


("C"  for  CEAREX)  . 


(54) 


For  u  »  o0,  it  behaves  as  l/o2,  and  for  o  <  w0,  it  remains 
approximately  constant.  Choosing  o0  corresponding  to  4X10'2 
cycles  per  hour  gives 

w0  »  2 nfa  =  2n(l.  11  xI0'7  Hz)  =  7.0x20'’  rad/ sec.  (55) 

Measurements  of  internal  wave  spectra  in  other  oceans  of 
the  world  at  lower  latitudes  have  indicated  a  characteristic 
1/f2  frequency  dependence.  Garrett  and  Munk  (1971)  formulated 
a  popular  empirical  model  based  on  these  observations  with 
such  a  frequency  dependence.  However,  the  AIWEX  spectra 
revealed  less  total  energy  in  the  arctic  internal  wave  field 
by  factors  of  15-30  than  in  fields  described  by  the  Garrett- 
Munk  model  (Levine,  et  al.,1987).  The  CEAREX  data  indicated 
a  wavefield  less  energetic  than  lower  latitude  fields,  but 
still  significantly  more  energetic  than  the  AIWEX  wave  field. 
The  CEAREX  data,  due  to  loration,  were  subject  to  large 
diurnal  tidal  motions  which  pumped  energy  into  the  internal 
wave  field.  In  addition,  under  ice  topography  (ice  keels)  and 
local  sea  floor  geography  (  i.e.  Yermak  plateau)  affected  the 
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internal  wave  field  (Stanton) .  In  contrast,  the  AIWEX  data 
were  taken  in  the  Canada  Basin  of  the  Arctic  Ocean  (Padman,  et 
al. ,1990)  and  experienced  far  more  quiescent  conditions. 

Arctic  internal  wave  fields  possess  different 
characteristics  from  lower  latitude  fields.  Some  suggested 
causes  in  the  Arctic  are  less  forcing  due  to  wind,  a  damping 
effect  caused  by  the  ice  cover,  lack  of  surface  waves,  and 
weak  large  scale  circulation  (Levine,  et  al.,1987).  Choice  of 
frequency  dependence  in  the  internal  wave  field  model 
critically  depends  on  the  location  that  is  to  be  modeled.  An 
arctic  model  will  use  a  1/f  dependence  for  the  velocity  power 
spectrum,  whereas  a  lower  latitude  model  will  use  1/f2.  The 
CEAREX  location  represents  a  special  situation  where  large 
tidal  forcing  causes  the  internal  wave  field  to  possess 
characteristics  similar  to  lower  latitude  fields. 

Garrett  and  Munk  (1971)  show  a  collection  of  data  on 
displacement  power  spectra  versus  wavenumber  obtained  from 
towed  measurements.  The  data  follow  a  1/k2  pattern  over  a 
full  four  decades.  The  largest  wavenumber  corresponds  to  a 
wavelength  of  10,000  m.  Assuming  a  simple  plane  wave 
description,  the  velocity  spectrum  will  possess  the  same  l/k2 
dependence 

(99 

displacement  ■  e1  <**-“'> ,  velocity  -  d displacement)  .  Jw  <L.#k 
Jk7  dc  Jk7 

A  simple  empirical  model  for  the  k  dependence  could  resemble 
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the  omega  dependence  model  with  a  power  limiting  k02  factor. 
However,  to  facilitate  later  mathematical  computations,  use 
instead 


The  total  energy  in  the  spectrum  is  now  limited  by  the  cutoff 
at  kmJn  instead  of  having  the  spectrum  level  out  at  low  k. 
Assume  this  model  applies  to  both  arctic  and  lower  latitude 
internal  wave  fields. 

Combining  both  forms,  a  model  for  the  velocity  spectra  for 
the  Arctic  (i.e.  AIWEX)  is 


12  = 


cl 


^g>20+u2  k2  2 njff+f2  k‘ 


W  k  <  °°*  (58) 


The  model  for  lower  latitude  and  CEAREX  velocity  spectra 
becomes 


12  = 


cl 


cl 


Kin<  k  <  «.(59) 


(q20  +  w2)  k2  4 n2(fl*f2)  k2 

Comparison  with  the  AIWEX  and  CEAREX  data  determine  the 
constants  of  proportionality,  CA  and  cc.  At  1  cph,  the  AIWEX 
spectrum  is  approximately  0.3  (cm/s)2/cph  (0.11  (m/s)2/Hz),  and 
all  the  CEAREX  spectra  have  approximately  the  value  0.5 
(cm/s)2/cph  (0.18  (m/s)2/Hz).  Approximating  the  AIWEX  model 
as  proportional  to  1/f  at  this  frequency 


27 


(60) 


Cl  r-  j_..  m  _Ca _ 1_ 

2nfJkniBk2  “  2n£  k^' 


Evaluating  at  1  cph 


|va(l  cph)  |2  =  0.3 


(cm/s) 2  s 
cph 


cl 


2n  (1  cph)  k^n ' 


(61) 


Then 


C\  =  2k  (0.3)  (cm/s)2  icmln  =  1.9x10'*  A„in  (JL) 

o  • 


(62) 


Approximating  the  CEAREX  model  as  proportional  to  1/f2  at  this 
frequency 


|vej2  =  - i— -  ~dk  ■  - i-—  — i— , 

c  4 n2fJ  ■»**,,  Jc2  4R2f2 


(63) 


Evaluating  at  1  cph 


Cc 


|vc(i  cph)  |2  =  0  5  (cm/s)2  m  — 
c  cph  4n2  (1  cph)2  lc^B 


(64) 


Then 

Cl  =  4n2  (0.5)  (cm/ s)2-cph  k min  =  5. 5x1  O'6  (  —  )  .  (65) 

s3 

The  formulas  for  the  induced  magnetic  field  involve  the 
velocity  amplitude  and  not  the  velocity  spectral  density.  The 
conversion  to  a  velocity  amplitude  was  made  using  the 
definition  of  the  power  spectral  density.  The  conversion  for 
this  case  was  simple  because  of  the  simple  time  dependence. 
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Given  a  function  f(t),  the  power  spectral  density  is  the 
Fourier  transform  of  the  autocorrelation  function  of  f(t). 
Let  F ( f )  be  the  Fourier  transform  of  f(t).  The  correlation 
theorem  says  (Brigham,  1988) 

autocorrelation  of  fit)  «  J  fi*)fit+t)dx  -  FT  pair  -  Fif)F'i/)**^ 

The  power  spectral  density  is  thus  F(f)F*(f).  Assuming  the 
velocity  has  the  following  time  dependence 

vit)  -  fit)  =  C  eiuc,  («7) 

then  the  Fourier  transform  of  the  velocity  is 

V(f)  =  r  C  ei2'nei2*t'tdt  =  C  S  if-f0)  .  (68) 


The  power  spectral  density  becomes 

\vif)\*  =  Vif)V’(f)  =  C2  5 2  if-f0)  =  C2 .  (69) 


Thus,  for  this  simple  time  dependence  and  real  coefficient  C, 
the  power  spectral  density  is  obtained  simply  by  squaring  the 
coefficient. 

Applying  this  result  to  the  above  model,  the  velocity 
becomes 


vA{k,u) 


(u.  ♦  u2) 


2  \  1/4 


v-ik,  w) 


fi?, 


Kin  <  *  <  - 


(70) 


♦ 
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This  simple  model  in  which  the  frequency  dependence  factor  and 
the  wavenumber  dependence  factors  are  multiplicative  factors 
makes  the  following  calculations  much  easier.  However,  as 
Garrett  and  Munk  (1971)  point  out,  actual  observations 
indicate  that  a  given  wavenumber's  contribution  depends  on  the 
frequency.  Following  Garrett  and  Munk,  the  function  m(u)  was 
included  here  as  a  reminder  of  this  additional  frequency 
dependence 


vA(k,u) 


CA  m.  (4l> ) 

<uJe  ♦  «2)1/4  k' 


vc(k,u) 


Cc  u  ( a) ) 

♦  o>*  k 


<  *  <  *• 


(71) 


Equation  [71]  replaces  the  constant  "C"  in  the  induced 
magnetic  field  expression,  Equation  [50]. 

Examine  two  different  methods  to  remove  the  k-dependence. 
For  the  first  method,  assume  internal  waves  exist  with  a 
continuum  of  frequencies  from  zero  to  the  buoyancy  frequency, 
N,  and  assume  they  possess  a  continuum  of  wavenumbers  from  a 
minimum,  kffljn,  to  infinity.  The  rationale  is  that  actual 
internal  wave  fields  result  from  a  complex  superposition  of 
waves  generated  by  different  sources  propagating  in  a  non- 
uniform  ocean.  The  assumption  of  a  continuum  may  circumvent 
a  detailed  description  of  this  wave  field.  The  second  method 
utilizes  the  eigenfunction  mode  structure  of  the  wave  field 
required  by  Equations  [17]  and  [19].  Application  of  a 
boundary  condition  on  the  vertical  velocity  at  depth  D  yields 
a  dispersion  relation  between  o  and  kx  for  each  mode.  These 
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dispersion  relations  remove  the  kx  dependence  in  B'.  The 
power  spectral  density  B'2  is  calculated  for  the  lowest  order 
mode  n=l .  The  lowest  order  modes  will  contain  a  majority  of 
the  energy  in  an  internal  wave  field. 

Try  the  continuum  approach  first.  To  simplify 
computations,  assume  that  the  depth  D  »  l/kx.  Then  the  last 
two  terms  in  the  induced  magnetic  field  expression,  which  are 
weighted  by  exp(-kxD)  ,  can  be  neglected  with  respect  to  1, 
simplifying  the  conversion  to  power  spectral  densities.  The 
induced  magnetic  field  has  the  simple  e’wt  time  dependence 
referred  to  above.  Squaring  the  amplitude  gives  the 
corresponding  power  spectral  density 


W4i/w20-*-w2 


(72) 


|Bc|2  =  [n  0BoCc)  2 


i2(ti))  fa)2  (N2  +  <j>2) 
N 4  (w20+g>2) 


The  formulas  above  are  "one-dimensional"  formulas  in  that 
they  were  derived  assuming  only  one-dimensional  wavevectors, 
kx.  The  more  general  wavefield  will  be  two-dimensional  with 
wavevectors  of  magnitude 

k .  ji&k*.  <™> 


The  one-dimensional  expressions  for  the  induced  magnetic  field 
power  spectrum  must  be  converted  to  a  two-dimensional  form. 
Yaglom  (1987)  shows  a  way  to  accomplish  this.  Given  a  one- 
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dimensional  power  spectral  density  function,  f,(k,).  Then  the 
corresponding  two-dimensional  power  spectral  density  function 
is 


f2(k) 


_1  p  d  dk1 

*  Jt  dk‘  fit  -  k>  ' 


(75) 


The  one  dimensional  induced  magnetic  field  k  dependence  was 
l/kx4.  Appendix  D  contains  a  summary  of  the  details  in 
computing  a  two  dimensional  k  dependence.  Then 

'i  <«  a  T  T?  •  <7‘> 


Now,  perform  the  integration  over  k  (kmin  <  k  <  «) 

1JL  dk _ Ll  r  -  3  1 

\lD  4  ic5  16  ic4  '**ln  16  j^n  ' 


(77) 


The  induced  magnetic  field  two-dimensional  power  spectral 
density  becomes 


[H0o  BCa)2 

(|i  0oBCc]2 


3 

16  kj, iin 

3 

16  ktn in 


H2((0)  (N2-0i2)O2 
u20+(j>2 

U2(o)  (N2~u>2)  <i?2 
N4  (w2e+w2) 


(78) 


(79) 


Now  apply  the  second  method  using  the  modal  structure. 
Combining  Equations  [26], [27]  and  [28]  gives  the  vertical 
component  of  velocity 

w  =  (C1eibz  *  Cze'ibz)  ei(fc'x'w£),  (for  ky=0)  .  (80) 
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Let  C,  ■  -C2  *  C/2i.  Then 

w=Csin(6z)  ei(k"* '  wt) .  (81) 

As  a  boundary  condition,  require  that  the  vertical  velocity 
equal  zero  at  the  lower  boundary  of  the  internal  wave  field  at 
depth  D.  Then 

w\,mD  =  C  sin(fiZ»  eilk**  •  "rt  .  o,  (82) 


sin(8l>)  =0,  —  bD  = 


N 


N2-<j)2 

o2 


kj)  =  nn ,  n= ? ,  2 


(83) 


kx  = 


nn 
D  \ 


o' 


N2-i o2  ' 


U2  = 


klD2N2 


n2n2 


k\D2 


(84) 


Replace  kx  in  the  formula  for  B*  in  Equation  [50]  with  this 
dispersion  relation.  Taking  h=0  and  mode  number  n=l  gives 

T  r^~ 

B'  =  (i 0oBC  e'i“t  — (w2-g>2)  [  l  +  e  V  w’-w1  ]  ,  (85) 

n/72 


Converting  to  a  power  spectral  density 


Bn  =  [p  gQB)2C2~—  (W2-u2)2  [  1  +  e'^  N“-“2 
n2N* 


j2. 


(86) 


Model  the  velocity  spectra  with  only  a  frequency  dependence 
since  a  dispersion  relation  exists  relating  kx  to  w. 
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(87) 


C\ 


|vc|s  - 


cl 


o  0  ♦  u* 


«o  +  tt2 


Replace  C2  in  the  B'2  expression,  Equation  [86],  with  Equation 
[87] 


Is'*!1  *  I  1  *  e"'**'  )’, 

yta3 


»20+<,)a  If3  If 4 


(88) 


tie  I2  *  [  n  0o  S] 2  — - £L(W2-w2)2  [  1  +  ]2> 


fa)2o+C02  ^2If* 


(89) 


C.  UNDER  ICE  TURBULENCE  MAGNETIC  SIGNATURE 

1.  Geometry 

The  same  geometry  as  for  the  internal  wave  calculation  was 
used,  except  that  now  the  y-coordinate  was  important  (Fig.  5) . 
The  distance  R  was  given  by 

|^|  *  lx2  +  y2  +  (h+z) 2 ,  with  R  =  x£  +  y/  +  z£. 

For  calculational  purposes,  turbulence  was  assumed  to  be 
confined  to  the  region  D0<z<D. 

2.  Turbulence  magnetic  signature 

The  expressions  for  the  earth's  magnetic  field  and  the 
induced  current  density  are 
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Figure  5:  Geometry  for  Calculation  of  Turbulence  Induced 

Magnetic  Field 
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B  =  B  fc,  and  J  =  o(v  x  B)  . 


(91) 


Due  to  the  vertical  direction  of  the  magnetic  field  in  arctic 
regions,  the  current  density  has  only  x  and  y  components 

J  *  o(VyBl  -  V^Bf)  .  (92) 

Make  a  change  of  variables  for  the  vertical  direction 

z+h—zsdz—dz;  range  D0<z<D  —  D0+h  <  z  <  D+h.  (93) 


Then 

J  x  R  *  aB[-vxz£  -  vyzf  ♦  (Vyy+VfX)  fi\  .  (94) 

The  Biot-Savart  law  then  gives 

*  (x*  ♦  y*  *  z2)  T 

Integration  over  x,y,z  gives  the  components  of  the  total 
induced  magnetic  field. 


Ho0B  fy*  rt’D*h _ -va _ 

4*  L-Jy.  Jr*De.»  fiJTkTjui+U2 


eiUcrx-ut)  z 

(j<3+y2+z2) 3/,J 


dxdydz . 


(96) 


The  y-component  equation  is  similar. 


t*  0°B  !•£>•<! (  xeJkwX 

4  ft  t*D9*h  ^  ^  2 


The  z-component  is 

— — - ° - dxdydz,  (97 ) 

(x’^y’+z1) 3/3 


Performing  the  integrations  [see  Appendix  B)  yields 
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(98) 


B 


/  _  -li0 0Bvoe-^ 


2kx)jkl+klju30+u>2 


e  6  e 


k*h  [e~k,fio  -  e'**0]  , 


B 


/  _  -\i0oBv0e 


■lut 


2kyJki+ky  J<ji20+(ji3 


,-i«t  g'V*  [g'*y°0  _  ^'V5]  , 


(99) 


B/_  =  i\i  0QBv0e‘iut  f  e-*^(e-*«gc-e-*^>)  +  e -V»(e g-V*)  1  .  (100) 
2iJo20*u2  kjjkl+kl  kyjkl+ky 


Approximating  the  geometry  for  an  on-the-ice  measurement 

h  -  0;  De  -  0;  D  *•*>.  (101) 


B 


-iuc 


x(y)  ~ 


2\f<y 


,+u‘ 


kx'jkl*kl{y 


(102) 


/  ip  0oBv0e*iuc 


b2  = 


2v/o>20+u2  kx\jkl+kl  kyJkl+k 


]  •  (103) 


Letting  D  approach  infinity  greatly  simplifies  the  following 
calculations.  This  procedure  is  justified  due  to  the  expected 
magnitude  of  the  parameters 

kx  =  »  2*  ;  Z>-100-1000m;  <  0 . 002  <  1.(104) 
*  A,  100m 

These  equations  still  have  a  k-dependence.  To  get  the 
components  as  functions  of  only  o  and  t,  integrate  over  the  k- 
dependence 
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1 


(105) 


,  -p  aoBv0e 

°K  * 


-i«e 


2 


kjjkl+k. 


dk . 


b  -\L0oBv0e-lut 

2^«e+«2 


1  r  Jkl+kl in  +  £0 

*o  *»in 


(106) 


where  k^  is  the  smallest  wavenumber  (corresponding  to  the 
longest  wavelength)  present  in  the  x-direction.  A  limit  on 
wavelength  is  realistic  for  the  ocean,  and  required  to  keep 
the  integral  finite.  Recall  that  ke  keeps  the  total  power 
finite  as  k  approached  zero.  k0  thus  represents  the  point  in 
the  power  spectrum  at  which  the  characteristic  k'z  behavior 
begins.  k0  is  thus  a  reasonable  choice  for  Jcmln.  Set  k^  equal 
to  k0.  Then 


B/  .  -li0oBv0e-ifa>c  lnly/T*l) 
2 


(107) 


Similarly  for  the  y  and  z  components 


B,_  -n  (1oBvce'lut  ln(^i)  Qnd  B,_  ip  0QBvoe’iue  ln(^l)  (MB) 


For  the  z-component,  integration  was  done  only  over  the  k 
components  which  gave  rise  to  the  respective  term  via  the 
Biot-Savart  formula. 
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The  total  induced  magnetic  field  is  now 


B'  •  \B,\»JbI2x  +  B^y  «■  BI3s • 


\io0Bvoe-lut  ln(s/?+l)  _T 
2jo30+<i>2 


(109) 


The  horizontal  component  of  the  induced  magnetic  field  is 

JB'\+Bny  -  -  -  .  (HO) 


B'h 


\i0QBv0e-lue  ln(^l) 


Again,  the  power  spectrum  was  obtained  by  Fourier 
transform  of  the  autocorrelation  of  B',  which  is  equivalent  to 
squaring  the  amplitude  h6re 

3  ( |i  o0Bvo) 2  <ln(v/F+l))2 


B 


/  2 


2  (oV“2)  kl 


(111) 


Bn  _  (»o°Bv0)2  ( In  (y/2~+l) )  2 

*  2  (v>z0+u>2)  kl 

D.  MEASURED  MAGNETIC  SPECTRA 


(112) 


The  Arctic  Internal  Wave  Experiment  (AIWEX)  was  conducted 
in  1985  on  the  arctic  ice  at  74 'N,  145 *W.  In  April  of  that 
year,  P.  Czipott  and  W.  Podney  used  two  induction  coil 
magnetometers  to  measure  the  North-South  horizontal  component 
of  the  Earth's  magnetic  field.  Fig.  6  shows  the  results  from 
three  runs.  The  spectra  flattened  out  at  low  frequencies 
because  the  electronic  gain  declines  at  low  frequencies.  The 
region  in  the  10*2  -10‘3  Hz  decade  was  relatively  unaffected. 
The  spectra  here  showed  a  cnaracteristic  1/f2  dependence. 


39 


They  also  showed  an  approximately  hundred- fold  magnitude 
variation  between  the  three  runs.  Thus,  large  magnitude 
changes  can  occur  over  a  temporal  scale  on  the  order  of  days. 
Note  that  the  bottom  spectrum  basically  coincides  with  the 
dashed  line  spectrum  which  corresponds  to  measurements  made  on 
land  at  Poker  Flat,  Alaska. 
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tall'/M/ 


niwvicr  (oh) 


Power  apectra  for  magnetometer  output  from  Rune  100 
(cantor)#  103  (top)#  and  107  (bottom)  at  the  lea  camp# 
and  Run  124  (daahod  curve)  froa  Pokor  Plat.  Tho  acalo 
la  calibrated  by  tho  magnetometer  gain  at  0.1  Bx; 
figure  3  llluatrataa  the  dropoff  in  gain  at  lower 
frequency.  Error  bar  ahowa  the  90%  confidence 
interval. 


Figure  6:  AIWEX  Power  Spectra  for  magnetometer  (Czipott, 

Podney,  1985) 
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IV.  RE8ULTS 


A.  EVALUATION  OF  INTERNAL  WAVE  MAGNETIC  FIELD  MODEL 

Numerical  evaluation  of  the  formulas  from  the  model  used 
the  following  values/parameters 

\k0  =  4 it xlO'1  -J^}  o  »  4.0  -jjj;  B  =  25,000  nTi  (113) 

-I2-- 3. 3 2x1 0'3  tf-1.05 xlO'2-^;  (114) 

az  m 4  sec 

=  1C)S  -  ^lin  *  6.28 (115) 

Table  II  lists  sample  numerical  values  using  the  continuum  w- 
and  k-dependence  method.  Table  III  lists  sample  numerical 
values  using  the  modal  method.  For  the  continuum  method,  Fig. 
7  shows  the  AIWEX  model  spectrum,  the  CEAREX  model  spectrum, 
and  an  extrapolated  spectrum  for  the  AIWEX  magnetometer  data. 
Extrapolation  was  made  from  the  10'1  -  10'2  Hz  range  of  Run 
100.  Fig.  8  shows  the  same  spectra  using  the  modal  method  for 
the  AIWEX  and  CEAREX  models.  The  vertical  extent,  D,  of  the 
internal  wave  field  was  300  m  (D'Asaro,  Morehead,  1991) . 
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B'*2  (nTA2/Hz) 


Frequency  (Hz) 


Figure  7:  AIWEX  Model  Spectrum,  CEAREX  Model  Spectrum,  and 
Extrapolated  Spectrum  for  AIWEX  Magnetometer  Data  (Continuum 
Method) 
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B'^2  (nT  ^  2/Hz) 


Figure  8:  AIWEX  Model  Spectrum,  CEAREX  Model  Spectrum,  and 
Extrapolated  Spectrum  for  AIWEX  Magnetometer  Data  (Mode 
Structure  Method) 
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Table  II:  Internal  Wave  Model  Numerical  Values 
(Continuum  Method) 


frequency  (Hz) 

10'5 

10'4 

10'3 

1 . 6xl0‘3 


AIWEX  JDT-2/H2) 

1.3X103 

1.3X104 

8.2X104 

1.7X104 


B.1  _  CEAREX  _  f  nT2/Hz ) 

6.0X104 

6.0X104 

3.8X104 

4.9X103 


Table  III:  Internal  Wave  Model  Numerical  Values 
(Mode  Structure  Method) 


f£egy,gn,gy,_£MZ) 

IQ'5 

10'4 

10'3 

1.6X10'3 


B1  AIWEX  (nT2/Hz) 

2.7X106 

2.7x10s 

1.1X104 

1.2X102 


B1  CEAREX  fnT2/Hz) 

1.3X108 
1 . 3x10° 

5.2X103 

3.3X101 


1 


The  AIWEX  magnetometer  measurements  were  not  compensated 
for  ionospherically  generated  magnetic  fields.  To  separate 
the  internal  wave  induced  magnetic  fields  from  the  ionospheric 
fields,  Czipott  and  Podney  placed  a  tiltmeter  on  the  ice  as  an 
independent  measure  of  internal  wave  activity  at  the  AIWEX 
site.  They  also  obtained  magnetograms  of  ionospheric  activity 
during  the  experiment  from  several  measurement  stations  in 
Alaska.  Following  Czipott  and  Podney,  the  signals  from  these 
three  sources  are  compared  in  Table  IV. 
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Table  IV:  Signal  Characterizations 


RUM  « 

IONOSPHERE 

TILTMETER 

MAGNETOMETER 

(MAGNETOGRAM) 

100 

"quiet" 

"high" 

"active" 

■•03 

"active" 

"active" 

"high" 

5 

"active" 

"quiet" 

"quiet" 

107 

"quiet" 

"quiet" 

"quiet" 

As  a  first  order  analysis,  assume  the  ionosphere  had  only  two 
levels  of  activity,  "quiet"  and  "active",  measured  by  the 
magnetograms.  Assume  the  under  ice  seawater  motion  also  had 
only  two  levels  of  activity,  "high"/"active"  and  "quiet", 
measured  by  the  tiltmeter.  Then  the  above  four  runs  show  the 
resulting  magnetometer  signal  for  the  four  possible 
combinations  of  ionospheric  activity  and  seawater  motion. 
Runs  100  and  105  indicate  th^  ;  the  ionosphere  and  the  seawater 
motion  each  produced  a  measurable  magnetometer  signal. 
Together,  runs  100,  103,  and  105  indicate  that  the  magnitudes 
of  the  magnetometer  and  ionospheric  signals  were  of  the  same 
order.  Run  107  affirms  that  the  magnetometer  readings  result 
from  actual  ionospheric  and  seawater  motion  signals. 

Extrapolation  of  the  AIWEX  magnetometer  data  with  a  1/f2 
frequency  dependence  to  the  10'3  -  10'5  Hz  region  gave  power 
spectrum  values  between  102  and  106  nT2/Hz.  For  the  continuum 
method,  the  model  values  varied  one  to  two  orders  of  magnitude 
above  and  below  these  extrapolated  values.  The  l/kmtn 
dependence  made  these  continuum  models  sensitive  to  the  value 


46 


of  k^  chosen.  Choosing  Jlmax  *  lot  )  m  instead  of  10,000  m 
reduced  the  values  by  an  order  of  magnitude.  For  the  modal 
method  retaining  only  the  lowest  order  mode,  the  AIWEX  model 
predicted  values  one  to  two  orders  of  magnitude  higher  than 
the  extrapolated  data,  while  the  CEAREX  model  predicted  values 
about  two  orders  of  magnitude  above  the  extrapolated  data. 
The  mode  structure  model  values  depended  on  vertical  extent  of 
the  internal  wave  field  as  D2.  Taking  D  «  100  m  instead  of 
300  m  reduced  the  spectra  by  one  order  of  magnitude  and 
brought  them  within  one  order  of  magnitude  of  the  extrapolated 
data.  The  vertical  velocity  profile  of  Equation  [81]  has  a 
simple  sin(«Sz)  dependence.  Refining  this  profile  could 
significantly  reduce  the  energy  in  each  mode  since  B'2  is 
proportional  to  (velocity)2.  Overall,  the  models  reasonably 
predicted  the  magnitude  of  the  seawater  induced  magnetic 
fields. 

The  frequency  dependencies  of  the  continuum  models  varied 
significantly  from  the  1/f2  dependence  of  the  extrapolated 
data.  The  AIWEX  model  predicted  a  linear  increase  with 
frequency  and  then  a  rapid  falloff  as  the  buoyancy  frequency, 
N,  is  approached.  The  CEAREX  model  predicted  a  parabolic 
falloff  with  frequency  as  (N2-o2) .  Unfortunately,  the  gain 
falloff  of  the  AIWEX  data  over  internal  wave  frequencies  may 
have  masked  any  dropoff  near  the  buoyancy  frequency.  These 
continuum  models  predicted  too  much  energy  at  the  higher 
frequencies.  Apparently  the  assumption  of  frequency  and 
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wavenumber  continua  lead  to  erroneous  predictions.  The 
specific  relation  between  o  and  k,  denoted  by  m2(u)  in  the 
models,  roust  be  determined.  If  one  assumed  that  m(«)  goes  as 
1/ (number  of  wavenumbers) ,  and  that  number  of  wavenumbers  was 
proportional  to  frequency  (see,  for  example,  Garrett  and  Munk, 
1972,  for  a  similarly  modeled  dependence),  then  the  models 
would  more  closely  predict  the  u'2  dependence  of  the  data.  One 
could  compare  the  internal  wave  model  in  the  10'1  -  10*2  Hz 
range  because  internal  wave  oscillations  are  limited  to 
frequencies  below  the  buoyancy  frequency. 

The  mode  structure  method  gave  models  with  a  much  more 
realistic  frequency  dependence.  The  AIWEX  model  declined  as 
1/f  before  going  to  zero  near  the  buoyancy  frequency.  This 
dependence  paralleled  the  observed  AIWEX  velocity  spectra 
frequency  dependence.  The  CEAREX  model  fell  off  as  1/f2,  just 
as  the  CEAREX  velocity  spectra  did.  These  parallel  frequency 
dependences  between  induced  magnetic  field  and  generating 
velocities  satisfied  physical  intuition.  This  close 
correspondence  between  model  and  data  implied  that  a  proper 
description  of  internal  wave  induced  magnetic  fields  requires 
utilizing  the  internal  wave  modal  structure. 

B.  EVALUATION  OF  TURBULENCE  MAGNETIC  FIELD  MODEL 

Numerical  evaluation  of  the  formulas  use  the  following 
values/parameters 
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Mo  =  4 piXlO'1-^-;  o  =  4.0 

A2  ^ 


vl  =  5.7XJ0'4  kr 


m 1 


B  -  25, OOOnT;  (116) 
to  _  =  2 ltXlO'*Hz.  (117) 


°  Hz-S 4 

The  v02  value  was  determined  the  same  way  as  C2  for  the 
internal  wave  model,  assuming  that  the  CEAREX  data  1/f2 
dependence  extrapolates  to  the  lO^-lO*2  Hz  range.  Sample 
numerical  values  using  Equation  102  are  given  in  Table  V. 
Fig.  9  shows  a  plot  of  the  BH 1  power  spectral  density  (nT2/Hz) 
versus  frequency  for  *  100  m  and  also  the  extrapolated 
AIWEX  magnetometer  spectrum . 

Table  V:  Turbulence  Model  Numerical  Values 

(\*x“100  n> 


FREQ  (HZ) 

B‘H*(nT2/HZ) 

10*3 

4.4X10'6 

10*2 

4 . 4xl0‘® 

5X10*2 

1.8X10*’ 

10*1 

4.4X10*10 

The  model  spectra  exhibited  a  characteristic  1/w2 
dependence,  just  as  the  AIWEX  data  show  in  the  10'1  -  io*2  Hz 
range.  However,  the  model  values  were  seven  orders  of 
magnitude  less  than  the  AIWEX  data.  At  10'3  Hz  the  model 
predicted  magnitude  was  approaching  the  limits  of  the  most 
modern  field  detectors  at  10'6  (nT2/Hz)  .  The  smallness  of  the 
predicted  signal  derived  from  the  small  velocities  used  and 
the  maximum  characteristic  wavelength  of  100  m.  The  small 
velocities  on  the  order  of  10*1  -  10‘3  (cm2/s2-Hz)  were 
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experimentally  measured  during  AIWEX  (McPhee,  1989)  and  their 
power  spectrum  is  shown  in  Fig.  10.  The  assumption  of 
isotropic  turbulence  in  the  model  restricted  wavelength  values 
to  less  than  about  100  m  because  the  maximum  observed  depth  of 
the  upper  boundary  layer  is  about  100  m  (Stanton)  (note:  the 
"average"  depth  is  about  30-40  m)  .  Thus,  while  turbulent 
motion  in  the  upper  boundary  layer  of  the  ocean  could  produce 
magnetic  signals,  the  magnitude  remained  small  because  the 
volume  of  water  generating  the  signal  was  limited.  This 
implied  that  the  AIWEX  magnetometer  measurements  in  the  10*1  - 
10*5  band  were  due  to  ionospheric  signal  only. 
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Figure  9:  Turbulence  Model  and  Extrapolated 

Magnetometer  Spectra 


AIWEX 
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29.96 

33.71 

33.46 


Power  spectra  of  the  vertical  velocity  time  series  at  each  level.  The 
95  percent  confidence  interval  is  shown  at  lower  left.  The  spectra  are 
plotted  on  log-log  axes,  with  successive  spectra  displaced  one  decade  down 
for  clarity. 


Figure  10:  AIWEX  Vertical  Velocity  Power  Spectra 
(McPhee,  1989) 


V.  CONCLUSIONS 


The  magnetic  fields  induced  by  the  motions  of  seawater  in 
the  earth’s  magnetic  field  were  modeled  analytically.  Such 
analyses  required  mathematically  modeling  the  velocities  of 
the  internal  wave  field  and  internal  turbulence  field,  and 
then  applying  Ampere's  law  via  the  Biot-Savart  formulation. 
For  internal  waves,  the  values  predicted  by  the  model  for  the 
induced  magnetic  field  were  roughly  equal  in  order  of 
magnitude  to  extrapolated  magnetometer  data.  Treating 
internal  wave  frequency  and  wavenumber  as  continua  yielded 
unphysical  frequency  dependences.  Using  an  eigenmode 
structure  for  internal  waves  yielded  more  realistic  frequency 
dependences.  Thus,  proper  analysis  of  internal  wave  magnetic 
signatures  must  incorporate  the  internal  wave  modal  structure. 
Additional  magnetometer  measurements  should  certify  the  1/f2 
frequency  dependence  and  determine  the  behavior  near  the  local 
buoyancy  frequency.  The  turbulence  model  successfully 
coincided  with  magnetometer  measurements  in  frequency 
dependence.  The  predicted  frequency  dependence  matched  the 
l/fz  dependence  of  the  data.  However,  the  predicted 
magnitudes  were  several  orders  of  magnitude  below  the  AIWEX 
magnetometer  measurements  and  approached  the  limits  of  the 
most  sensitive  magnetometers.  This  implied  that  the 
relatively  small  scale  of  turbulence  motions  (lmax  =  100  m) 
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resulted  in  a  small  induced  magnetic  field  signature.  The 
AIWEX  data  in  the  10*1  -  10‘s  Hz  range  were  interpreted  as 
deriving  from  ionospherically  generated  signals  . 
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VI .  RECOMMENDATIONS 


Additional  measurements  of  seawater  motion  induced 
magnetic  fields  are  needed.  Such  measurements  should  be  made 
with  highly  sensitive  magnetometers,  for  example, 
superconducting  quantum  interference  devices  (SQUIDS) . 
Discrimination  against  ionospherically  induced  magnetic  fields 
must  be  done  by  using  the  large  spatial  coherence  length  of 
ionospheric  fields  to  identify  and  remove  their  signal. 

The  models  can  be  extended  to  predict  the  signal  seen  by 
a  gradiometer.  Thi^.  involves  expressing  the  induced  magnetic 
field  at  the  field  point  in  terms  of  field  point  position  and 
taking  the  gradient.  Gradient  measurements  of  internal  wave 
induced  magnetic  fields  are  more  numerous  than  magnetometer 
measurements  (for  example,  Czipott  and  Podney,  1985;  Podney, 
1975) . 

Both  the  internal  wave  and  the  turbulence  models  for 
induced  magnetic  field  utilized  the  approximately  vertical 
orientation  of  the  earth's  magnetic  field  in  the  Arctic  to 
simplify  the  calculations.  The  models  can  be  extended  to 
apply  the  Biot-Savart  calculation  method  to  the  non-vertical 
magnetic  field  at  other  parts  of  the  world's  oceans. 

The  internal  wave  model  assumed  a  constant  buoyancy 
frequency,  N,  for  the  entire  wave  field  in  the  horizontal  and 
vertical  dimensions.  In  the  actual  ocean,  the  buoyancy 
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frequency  changes  from  location  to  location  and  varies 
significantly  with  depth.  A  more  refined  model  must  include 
such  variations  and  should  predict  a  less  defined  cutoff 
frequency  due  to  the  spread  in  buoyancy  frequencies. 

The  modal  version  of  the  internal  wave  model  should  be 
refined  by  including  the  contributions  of  higher  order  modes. 
The  resulting  correction  should  not  be  too  large  because 
energy  content  rapidly  declines  for  higher  order  modes  due  to 
dissipation.  The  velocity  profiles  should  also  be  refined. 
As  an  example,  require  velocity  continuity  at  the  boundaries 
of  the  internal  wave  field  and  allow  exponential  decay  outside 
these  boundaries,  analogous  to  going  from  an  infinite  square 
well  solution  to  a  finite  square  well  solution.  The  amplitude 
within  the  well  is  reduced. 
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VII.  APPENDICES 


A.  CALCULATION  07  INTERNAL  WAVE  INDUCED  MAGNETIC  FIELD 

For  a  region  with  constant  Brunt-Vaisala  frequency,  N,  the 
vertical  component  of  the  seawater  velocity  is  given  by 

w(x,y,z,t)  =  1 7(z)  .  (118) 

p(z)  *  e~ik,z  [qe"*  ♦  C2e-i4']  .  (n») 


where: 


b 


N 


N2-l>2 

G>2 


(120) 


Require  as  a  boundary  condition 


w(x,y, z-o,  O «o,  —  (z»o)  *1  [q  +  c2] ,  —  ca«=-c2«c. 


(121) 


Then  the  vertical  velocity  component  becomes 

i-(z)  =  Cfc"1M[e"*  -  e'ii2]  .  d22> 

The  derivative  with  respect  to  z  is 

-4- it-)  -  C{- ik  +  ib)  e{~ik,tii) *  -  C(-ik.-ib)  e  (123) 

dz  1  * 

With  the  simplified  geometry,  the  y-component  of  velocity  is 
zero.  The  x-component  of  velocity  becomes 
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(124) 


u(x,y, z,  t) 


gi (lt-r  -  we) 


l£L  [©i«*+e'i4*]  ei{k**~*  } . 

K 


(125) 


These  expressions  can  now  be  used  to  calculate  the  induced 
magnetic  field 


B 


/ 

X 


Iff 


(126) 


J [•*»♦•  r  fD 
x»-“  Jz*D0 


Ho  Q(vBk  -  uB,) 
271  tf2 


(h+z) 


dxdz 


(127) 


f-rt-c  Bn0oC 
2  ft 


(s14*-«‘/4*)  cosycoa*  »  AsiM  (128) 

( - — -j - P - [*♦*)  *zdx 

<h»2)’  ♦  XJ 


The  x-integration  is  done  via  contour  integration  in  the  upper 
half  plane 


eu«* 

(h*z) 1  *  z 2 


dx 


l 

J  c 


W  (x*i  {A«-z) )  (x-i  (A»z) ) 


dx 


,JM 


(A+z) 3*zJ 


dx 


(129) 


=  27ti  [  — 


,  ikx(l(h*z) ) 


i  (h+z)  +  i  (h+z) 


] 


7t  n-kx(h+z)  (130) 
'  • —  S  • 

(h+z) 


Then 
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a 


r f * t •  •-« -) cosycos*  ♦  ISiBl  (,<•».»<«> S0 ».».-<»» »)dt 

Jo,  kM 


(131) 


Performing  the  z-intearations  and  noting  limits  of  integration 
yields 


B' 


- ?fq<~  e'*x/le'iMt  t  (  (-t—  r > e(~k*'ih)t+  (— -1_)  e <**«*•»>«) cosycos^  + 
2  Jo 


032) 


♦  jt ) > )  e 1  •*«-*♦>«]  j  |®.  (133) 

The  z-component  of  the  induced  magnetic  field  is 
calculated  similarly 


B‘ .  r**rc  Mx  .  f~r 

Jjcv"  Jt»P9  2 71  J?*  Jf)# 


V®  0  (^Bx  -  uBt) 


2  K 


x  dzdx 


(134) 


J— Jp.  2n  jeM  (A»z)J  ♦  x1 

«  1Hz?.aC 0-'-c[Dl  ( «««-«*«»)  co8YCos4*.££i2£  (#i*.+e-i*»)]  <i«*-**(A-*>)cte(136) 

275  »P#  Ajj 

«:  •  l  I  ( ( )  e '•*.•«>»)  eo6Tcoe»  *(137) 

*  l-jrht'  )  IS.-  (138) 


Note  that  the  x  and  z  components  have  the  same  magnitude  and 
differ  in  phase  by  ir/2.  These  last  equations  can  be 
simplified  by  multiplying  the  numerator  and  denominator  by  the 
complex  conjugate,  and  using  Euler's  identity  for  complex 
exponentials,  resulting  in 

iax-Bx=  ^  °^0C  ( — ~ — )  e-Jue  [  (-/c^inJz  -  Acos&z;  2icosycos$ 

2  kl+b2 


+  <-/cxcos6z  +  fisin  bz)  3  IS.-  (140) 
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B.  CALCULATION  OF  TURBULENCE  INDUCED  MAGNETIC  FIELD 

Tha  integrations  for  the  components  of  the  turbulence 
induced  magnetic  field  are  shown  below. 


»'  =  Ho oBfx-ry.-fM.o.j,  _-v0e'ivt _ g  ^  dxdydz  (141) 

4K  L-Jy - LDo.h  P^T/Wa+U2  (X2+y2  +  Z2)3/i 


'oTftA;V“  9’ 
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Jz-D.O  Jx—  {x+iz)  (x-iz) 


/“x«»  /*z«i>/>  * 


~2^  r,/0e' 
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Jx‘->Jz‘D0'h  (x2  +  Z2) 


Perform  x-integration  in  upper  half  of  complex  plane 

,  _  -2n.»av.  .  /•*■»•* 
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7i  f  e  ^  dz, 
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(145) 
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2kJkl  +  kl 


(146) 


The  integration  for  the  y-component  is  performed  in  the  exact 
same  manner,  with  x  and  y  interchanged. 

The  z-component  integration  is  calculated  similarly 

g,  mjr!cB^yr..r..^l _ V^xe"-* _ _ 1,(147) 

in  !JTv7T7kJ(x,*yl‘Z>),/J 
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c. 


INTERNAL  WAVE  VELOCITY  POWER  SPECTRAL  DEN8ITY  MEASURE  KENT  8 
Data  from  the  Coordinated  Eastern  Arctic  Experiment 


(CEAREX) ,  March-April  1989. 


Rotary  Spectro  50m  ADCP 


cm/s)Vcph  (cm/s)*/cph 


104 


Rotary  Spectro  100m  S4 


frequency  (cph) 
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Rotory  Spectra  150m  S4 


Rotary  Spectro  200m  S4 


frequency  (cph) 
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D.  I'D  TO  2-D  POWER  SPECTRUM  CONVERSION  CALCULATION 


The  one-dimensional  power  spectrum  and  its  derivative  are 


d  fjdcj) 
dkx 


(153) 


The  two-dimensional  power  spectrum  is  given  by 
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Make  a  change  of  variable 
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From  CRC  Standard  Math  Tables  (1987),  given  a  trinomial  in  m, 
M  =  a  +  bm  +  cm2,  then 


f  Mldm  =  -i-m”-W  - 
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Here, 

a=l,  b= 

f2(k) 
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Applying  this  integral  twice  yields 

Z  -  iZLHI  .  2sin‘*m  IS,  U59) 
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(160) 


67 


LIST  OF  REFERENCES 


Beal,  H.  T.  and  Weaver,  J.  T.,  “Calculations  of  Magnetic 
Variations  Induced  by  Internal  Ocean  Waves,"  Journal  of  Geophysical 
Research,  v.  75,  no.  33,  pp.  6486-6852,  20  November  1970. 

Berg,  P. ,  and  McGregor,  J.,  Elementary  Partial  differential  Equations , 
Holden-Day,  1966. 

Boas,  M.  L. ,  Mathematical  Methods  in  the  Physical  Sciences ,  2nd  ed.,  John 
Wiley  &  Sons,  Inc.,  1983. 

Brigham,  E.  o..  The  Fast  Fourier  Transform  and  its  Applications ,  Prentice- 
Hall,  Inc.,  1988. 

Chave,  A.  D. ,  "On  the  Electromagnetic  Fields  Induced  by 
Oceanic  Internal  Waves,"  Journal  of  Geophysical  Research ,  v.  89,  no. 
C6 ,  pp.  10519-10528,  20  November  1984. 

Condon,  E.  ,  and  Odishaw,  H.  ,  Handbook  of  Physics ,  McGraw-Hill, 
1967. 

CRC  Handbook  of  Chemistry  and  Physics,  62nd  ed.,  pp.  F-3,F- 
166,  CRC  Press,  Inc.,  1981. 

CRC  Standard  Mathematical  Tables ,  28  th  ed. ,  CRC  Press,  Inc.,  1987. 

Crews,  A.  and  Futterman,  J.,  "Geomagnetic  Micropulsations  Due 
to  the  Motion  of  Ocean  Waves,"  Journal  of  Geophysical  Research ,  v.  67, 
no.  1,  pp.  299-306,  January,  1962. 

Czipott,  P.  V.  and  Podney,  W.  N.,  "Measurement  of  magnetic 
Fluctuations  above  the  Arctic  Ocean  during  Project  AIWEX: 
Experiment  Description  and  Data  Overview,"  technical  report 
for  Physical  Dynamics,  Inc.,  1985. 

D'Asaro,  E.,  and  Morehead,  Mi,  "Internal  Waves  and  Velocity  Fine  Structure 
in  the  Arctic  Ocean" ,  Journal  of  Geophysical  Research,  v.  96,  no. 
C  7  ,  pp.  1  2  ,  7  2  5  -  12  ,  7  3  8  ,  July  15,  1991. 

Dietrich,  Kalle,  Krauss,  and  Siedler,  General  Oceanography ,  2nd 
ed.,  p.  62,  plates  3-5,  John  Wiley  &  Sons,  Inc.,  1980. 


68 


Garrett,  c.  and  Munk,  W. ,  "Space-Time  Scales  of  Internal 
Waves,"  Geophysical  Fluid  Dynamics ,  v.  2,  pp.  225-264,  1972. 

Gill,  A.  E.,  Atmosphere-Ocean  Dynamics,  Academic  Press,  Inc.,  1982. 

Gradshteyn,  I.  S.  and  Ryzhik,  I.  M. ,  Table  of  Integrals,  Series,  and 
Products,  Academic  Press,  Inc.,  1965. 

Gray,  D.  E.,  American  Institute  of  Physics  Handbook ,  McGraw-Hill,  1972. 

Gross,  M.  C.,  Oceanography,  A  View  of  the  Earth ,  pp.  156-160,  Prentice- 
Hall,  Inc.,  1972. 

Kraichman,  M.  B. ,  Handbook  of  Electromagnetic  Propagation  in  Conducting 
Media,  U.S.  Government  Printing  Office,  Washington,  DC,  1976. 

Levine,  M.,  Paulson,  C. ,  and  Morison,  J.,  "Observations  of  Internal 
Gravity  Waves  Under  the  Arctic  Pack  Ice ",  Journal  of  Geophysical 
Research,  v.  92,  no.  Cl,  pp.  779-782,  January  15,  1987. 

McPhee  M. ,  " A  Rigid,  Cable-Lowered  Instrument  Frame  for  Measuring  Turbulence 
and  Internal  Waves  in  the  Arctic" ,  IEEE  Journal  of  Oceanic  Engineering, 
V.  14,  no.  2,  pp.  203-207,  April  1989. 

Office  of  Naval  Research,  Data  Report  152,  Observations  from 
CEAREX  "O”  Camp,  Arctic  Ocean,  March-April  1989 ,  by  Levine,  M. , 
Paulson,  C.,  Simpkins,  J. ,  and  Gard,  S.,  Reference  91-1,  April 
1991. 

Parkinson,  Introduction  to  Geomagnetism,  pp.  72-76,  University 
Press,  1983. 

Peerless,  S.  J.  ,  Basic  Fluid  Mechanics ,  Pergamon  Press,  1967. 

Petersen,  R.  A.  and  Poehls,  K.  A.,  "Model  Spectrum  of  magnetic 
Induction  Caused  by  Ambient  Internal  Waves,"  Journal  of  Geophysical 
Research,  v.  87,  no.  Cl,  pp.  433-440,  20  January  1982. 

Phillips,  O.  M. ,  The  Dynamics  of  the  Upper  Ocean,  Cambridge 
University  Press,  1969. 

Podney,  W.,  "Electromagnetic  Fields  Generated  by  Ocean  Waves," 
Journal  of  Geophysical  Research ,  v.  80,  no.  21,  pp.  2977-2990,  20  July 
1975. 


69 


Podney,  W.  and  Sager,  R. ,  "Measurement  of  fluctuating  magnetic 
gradients  Originating  from  Oceanic  Internal  Waves,"  Science,  v. 
205,  pp.  1381-1382,  28  September  1979. 

Sommer f eld,  A.  W. ,  Electrodynamics,  Lectures  on  Theoretical  Physics,  v.  2, 
pp.  280-283,  Academic  Press,  Inc.,  1952. 

Stanton,  T. ,  Adjunct  Resident  Professor,  Oceanography,  Naval 
Postgraduate  School;  conversation  with,  11  December  1991. 

Tatar  ski,  v.  I.,  Wave  Propagation  in  a  Turbulent  Medium ,  McGraw-Hill 
Book  Company,  Inc.,  1961. 

Warburton,  F.  and  Caminiti,  R. ,  "The  Induced  Magnetic  field  of 
Sea  Waves ,"  Journal  of  Geophysical  Research ,  v.  69,  no.  20,  pp.  4311- 
4318,  15  October  1964. 

Weaver,  J.  T.,  "Magnetic  Variations  Associated  with  Ocean 
Waves  and  Swell,"  Journal  of  Geophysical  Research ,  v.  70,  pp.  1921- 
1929,  1965. 

Yaglcm,  A.  M. ,  Correlation  Theory  of  Staiionaiy  and  Related  Random  Functions, 
Vol.  It  Basic  Results ,  Springer-Verlag,  1987. 


70 


INITIAL  DISTRIBUTION  LIST 


1.  Library,  Code  52 

Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

2.  Prof.  Donald  L.  Walters 
Department  of  Physics  (Code  PH/We) 

Naval  Postgraduate  School 
Monterey,  CA  93943-5004 

3.  Prof.  Xavier  K.  Maruyama 
Department  of  Physics  (Code  PH/Mx) 

Naval  Postgraduate  School 
Monterey,  CA  93943-5004 

4.  Prof.  Karlheinz  E.  Woehler 
Department  of  Physics  (Code  PH/Wh) 

Naval  Postgraduate  School 
Monterey,  CA  93943-5004 

5.  Prof.  Timothy  P.  Stanton 
Department  of  Oceanography  (Code  OC/St) 
Naval  Postgraduate  School 

Monterey,  CA  93943-5004 

6.  Dr.  Harvey  W.  Ko 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  Md  20723-6099 


2 


5 


1 


1 


1 


1 


71 


